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Abstract 

We present a first-principles method for the calculation of optical excitations in nanosystems. 
The method is based on solving the Bethe-Salpeter equation (BSE) for neutral excitations. The 
electron self-energy is evaluated within the GW approximation, with dynamical screening effects 
described within time-dependent density-functional theory in the adiabatic, local approximation. 
This method is applied to two systems: the benzene molecule, C§R§, and azobenzene, C12H10N2. 
We give a description of the photoisomerization process of azobenzene after an n — ir* excitation, 
which is consistent with multi-configuration calculations. 
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I. INTRODUCTION 



Predicting many-body excitations of weakly correlated systems from first principles is 
important, both to help understand the behavior of known systems as well as to design 
novel materials and technological devices. The study of optical excitations in systems at 
nanoscale size has been accentuated by recent achievements in nanotechnology and the 
potential growth of the field . On the theory side, time-dependent density functional 

theory in the adiabatic, local approximation (TDLDA) has been used to study a wide class 
of nanosystems ranging from organic molecules to semiconducting and metallic clusters (see 
e.g. j^J and references therein). Alternatively, optical excitations of atoms, small molecules 
and bulk materials have been studied from first principles by solving the Bethe-Salpeter 
equation (BSE) for electrons and holes lal- In bulk materials, the BSE approach is known 

. n n 

to fully describe excitonic effects, which are missing in TDLDA [5|, |6J. In confined systems 
such as clusters and isolated molecules, the BSE approach is expected to be more accurate 
than TDLDA, but extensive comparisons of both methods have not been done so far due 
to the complexity of most numerical implementations of the BSE. One source of complexity 
is the explicit evaluation of the dielectric function of the system. Early implementations 
use a Fourier expansion of the static dielectric function, with dynamical effects included 
via semi-empirical models or ignored altogether. More recent studies employ generalized 
olasmon-pole models in the description of dynamical screening Ogiit and collaborators 
2| have presented a fully ah initio technique for the calculation of the static dielectric 
function in real space, leading to important simplifications compared to Fourier-expansion 
techniques. 

We propose a formulation of the BSE method that avoids the explicit evaluation of the 
dielectric function. All the information about electronic screening is contained in the po- 
larizabilitv operator, which is calculated within TDLDA by solving a generalized eigenvalue 



probtem flu After — dia^on, ei„ S g ,ve tfte spatial defence 
of the polarizability, while eigenvalues give its frequency (time) dependence. Both the elec- 
tronic self-energy and electron-hole interaction kernel, which enter explicitly in the BSE, are 
computed directly from the TDLDA polarizability. Once the BSE is reduced to a general- 
ized, frequency-dependent eigenvalue problem, it is numerically diagonalized. The resulting 
normal modes are associated to neutral, optical excitations of the electronic system. 
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II. THEORETICAL METHOD 



The underlying description of the electronic system is obtained within density-functional 
theory (DFT). The Kohn-Sham equations are solved in real space. Degrees of freedom 
related to core electrons are excluded from the problem by using norm-conserving pseudo- 
potentials. We start by expressing the random-phase approximation (RPA) polar izability 
operator in the space of Kohn-Sham eigenstates in the usual way: 



uo — oj cv + i0 + to + uo cv — i0 + 
where u cv is the difference between two Kohn-Sham eigenvalues, e c — e v , and + denotes an 
arbitrarily small quantity. We assume h = 1 and ignore spin indices for the sake of clarity. 
The case of spin-polarized systems require slight modifications in the formalism. Here, we 
denote occupied Kohn-Sham orbitals by v, v' and unoccupied (virtual) ones by c, c'. Generic 
orbitals are denoted by letters i,j,n. Appropriate occupation factors should be included in 
Eq. (^) if partially populated orbitals are present jjj. Within TDLDA, the polarizability 
is written in terms of eigenvalues {u m } and eigenvectors {F™} of an effective eigenvalue 
problem 0, 



^cv,c'v' / / 
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(2) 

LU — LU m + i0 + LU + LU m — iO -1 " 

Both polarizabilities can also be expressed in real space by including Kohn-Sham eigenvectors 
according to the rule 

{O} (r, r') = £ u,(rK(r') K>h,„ «$M«;'(*0 ■ (3) 

ij,i'j' 

One key quantity in solving the Bethe-Salpeter equation for optical excitations is the 
electronic self-energy E, which can be computed from first principles within the so-called GW 
approximation (GWA) This method has been used in a number of different systems 

and at different levels of sophistication. At its lowest level, the self-energy is schematically 
evaluated as S = G W , where G is the (DFT) Green's function and W = [1 — VPqI^V is 
the screened Coulomb interaction, with V being the bare Coulomb interaction, V(r) = e 2 jr. 
In this implementation, we go beyond RPA and use the TDLDA polarizability, Eq. (E). 
Having in mind that n and Pq are related to each other via n" 1 = Pq 1 — (V + f) p, Il2| . 
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where / is the TDLDA interaction kernel, / = ^r^, one can arrive at this expression for the 
self-energy: 



£(u/) = i / 0+ G (uj' - u) 



V + VU(u)V + \vil{u)f + \fTl{u)V 



(4) 



In Eq. (j4j), the first term inside square brackets denotes the bare exchange (Hartree) 
contribution. The second term is the correlation part. The last two terms have vertex 
corrections. Some aspects of Eq. (JH) should be pointed out: 

1. The additional vertex terms are written in symmetrized form, so that the resulting 
self-energy is symmetric with the interchange of arguments. In many-body notation, 
£(1,2) = £(2,1). 

2. By ignoring the / kernel and evaluating the polarizability as IT -1 = Pq 1 — V, we recover 
the self-energy at £ level of approximation. Adding a / kernel in the polarizability 
amounts to enhanced screening, which is partially compensated by the inclusion of 
vertex terms. 

3. The frequency dependence of the polarizability is known from Eq. (J2J). Plasmon pole 
models are not needed. 

Numerically, matrix elements of the self-energy can be evaluated by expressing the quan- 
tities V, II, / in representation of Kohn-Sham orbitals using Eq. (j3J) and evaluating the 
self-energy directly from Eq. (J3J). We do not impose self-consistency between polarizability, 
self-energy and Green's function. 

Optical excitations can be obtained either from the singularities of the TDLDA polariz- 
ability in frequency domain, Eq. (J2J), or by solving the BSE. In the second case, we set up 
a generalized eigenvalue problem where the interaction kernel K is no longer of the form 
K = V + f as in TDLDA, but given by K = ^+ s) Q. From Eq. ©, this reduces to 



Kvc,v'c' — Wlvt-v'cf 



VII f flLV 

v + vnv + + J 



(5) 
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With this kernel operator, the BSE is solved numerically using standard procedures. The 
erm inside the first pair of square brackets should be ignored for spin-triplet excitations 
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Having a frequency- dependent kernel makes the BSE a non-standard eigenvalue equation. 
In periodic systems, dynamical effects were found to be negligible if one is interested in the 
linear optical spectrum. In finite systems, quantum confinement effects increase the strength 
of the interaction kernel, resulting in somewhat stronger but still limited dynamical effects 
(about 0.01 eV for the SiH 4 system) Although we ignore those effects in the present 
work, they can be included by evaluating polarizability and kernel at some frequency close 
to the actual excitation frequencies obtained by solving the BSE. 



III. APPLICATIONS 



A test case of the procedure above is the isolated benzene (CqHq) molecule. For this 
system, we have solved the Kohn-Sham equations on a regular grid, with grid spacing 0.4 
a.u. (0.21 A). Electronic wave-functions were required to vanish outside a sphere of radius 
16 a.u. (8.47 A), centered on the molecule. Carbon-carbon and carbon-hydrogen bond 
lengths were fixed at their experimental values. After including self-energy corrections, 
the ionization potential is found to be 9.30 eV, in good agreement with the experimental 
value of 9.3 eV Q]. Benzene has a negative electron affinity of -1.12 eV, measured in 
resonance scattering experiments Q], which is in excellent agreement with the calculated 
electron affinity of -1.11 eV. Here, the negative electron affinity is interpreted as energy 
of the resonant state of the anion, and not the lowest energy required to extract one 
electron from the anion (which is zero because the anion is unstable). We observed that 
the inclusion of vertex corrections together with a TDLDA polarizability is essential for an 
accurate prediction of the ionization potential and electron affinity of benzene. In fact, the 
S approximation predicts those quantities to be 9.84 eV and -0.54 eV respectively. Defining 
a "HOMO-LUMO gap" as the difference between ionization potential and electron affinity, 
the present approach and the S approach agree in the value of the gap: 10.4 eV. This is 
somewhat consistent with the observation by del Sole and collaborators [l5j], who conducted 
a similar analysis in bulk silicon and found significant shifts in the valence band maximum 
and conduction band minimum, but only small change in the energy gap. 

The low-energy optical spectrum of this molecule is dominated by a ir — -k* complex, which 
arises from transitions between the degenerate highest occupied molecular orbital (HOMO) 
and lowest unoccupied molecular orbital (LUMO). Table H] shows a comparison between 
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TDLDA and BSE predictions for some excitations in the n — 7r*-complex. Although singlet 
transitions B\ u and E\ u , which are the dominant ones in the low-energy part of the absorption 
spectrum, are equally well described by both methods, there is a small but well-defined 
blue shift of dark transitions Bf u and B\ u within TDLDA. With the exception of the E\ u 
transition, excitation energies obtained by solving the BSE are typically underestimated with 
respected to measured quantities. The overall deviation between measured and calculated 
transition energies is 0.1 to 0.3 eV. 

The absorption spectrum of benzene, Fig. shows the bright E\ u transition, at 7.0 eV. 
The low, flat feature in the 6.0-7.0 eV range is due to coupling between n — it* transitions 
and vibrational modes 13, ^||, and it is absent in the calculated spectra because of the 
assumed structural rigidity. Beyond 10 eV, a number of sharp features in the measured 
spectrum results from transitions involving Rydberg states jl^i ll t| . The limited numerical 
accuracy in that energy range prevents a detailed identification of such transitions in the 
calculated spectrum. 

A more interesting system is the azobenzene (C12H10N2) molecule. This molecule has 
attracted considerable attention recently due to exciting experiments of light-driven me- 
chanical manipulation at nanoscale. Some of these experiments are: demonstration of light- 
driven manipulation of liquid droplets on an optically active surface jl|; energy conversion 
in a polymeric azobenzene chain 3]; synthesis of an azobenzene compound which exhibits 
hinge-like motion when photo- irradiated 3|. The basic phenomenon explored in these ex- 
periments is photo-isomerization: the structure of azobenzene is induced to change from a 
trans-azobenzene configuration (TAB) to a metastable cis-azobenzene configuration (CAB) 
upon absorption/emission of radiation with the appropriate frequency. On the theory side, 
the detailed structural distortion involved in the process has been subject of discussion. 
Early analyses based on a restricted configuration interaction (CI) calculation indicated 
that photo-isomerization following excitation to the first excited state occurs via inversion of 
one nitrogen in the azo group, causing an increase of 120° in the C-C-N bond angle ("inver- 
sion path"). On the other hand, recent multi-configuration and TDLDA calculations [lflll^ 
indicate that the most favorable path is a torsion of the N-N bond around its axis, causing a 
180° rotation of one phenyl group with respect to the other ("rotation path"). Diau has pro- 
)osed an alternative "synchronous inversion path" that can also lead to photo-isomerization 

3- 
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We have studied optical excitations in the isolated azobenzene molecule, in order to 
test the accuracy of the present BSE approach. Isomerization along the rotation path was 
simulated by initially fixing the CNNC dihedral angle and minimizing the DFT total energy 
of the system with respect to all other degrees of freedom. A sequence of configurations 
was then obtained for various choices of the dihedral angle ranging from 0° (CAB) to 180° 
(TAB), as discussed in j^oj. The Kohn-Sham equations were solved on real space, using a 



0.21 A). At the equilibrium 
I, which is consistent with 



boundary radius of 20 a.u. (10.58 A) and grid spacing 0.4 a.u. 
configuration, trans-azobenzene has ionization potential 8.5 eV 
the corresponding values of 8.6 eV, obtained from GW with a self-energy given by Eq. (j3J). 

The calculated potential energy of the system in its ground state and the two lowest 
excited states is shown in Fig. |21 At each configuration, the potential energy is computed 
as a sum of ground-state potential energy and excitation energy, the latter computed either 
within TDLDA or BSE. Calculated excitation energies for the CAB and TAB configurations 
are shown in Table ITT1 Deviations between the BSE and experimental excitation energies are 
small, within 0.3 eV. The CAB titt* excitation energy has larger deviation, of almost 0.8 eV, 
which may indicate a difficulty in the theoretical description of higher-energy excitations. 
On the other hand, excitation energies obtained within TDLDA are systematically lower 
than the ones obtained within BSE by more than 0.4 eV, and substantially different from 
measured excitation energies. 

Figure El shows that, once the system absorbs radiation and is promoted to its first 
excited state at either CAB or TAB configurations, there is no energy barrier preventing 
fast relaxation to an intermediate configuration, with dihedral angle around 90°. From that 
configuration, the system can decay to ground state and relax back to either the CAB or 
TAB configurations. This mechanism allows trans — ► cis or cis — ► trans isomerization 
after excitation to state ror* (the first excited state). Excitation to the second excited state 
is not found to induce isomerization, unless there is decay to state rwr*, since the profile of 
potential energy is rather flat along the rotation path. 

Photo-isomerization involving higher-order excited states, which has been observed in lab- 
oratory, can also be studied once four-particle and more complex excitations are included 
in the theory. Multi-configuration calculations [rol . [^J indicate that an excited state com- 
posed of two electrons promoted from the HOMO to the LUMO has a barrier-less potential 
energy profile along the rotation path. That excited state is not accurately described within 
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either TDLDA or BSE approaches in the current two-particle approximation. Despite this 
limitation, excited states of azobenzene can be studied within the BSE approach or, to a 
lower level of accuracy, TDLDA. 



IV. CONCLUSION 



In summary, we have discussed a first-principles method for the calculation of opti- 
cal excitations in nanosystems, with applications to two organic molecules: benzene and 
azobenzene. The method is based on solving the Bethe-Salpeter equation (BSE) for neu- 
tral excitations, with dynamical screening described within TDLDA. To our knowledge, this 
is the first time that optical excitations in organic molecules are calculated using Green's 
function-based methods. The approach is completely system-independent and it can also be 
used to study inorganic systems. This work was supported in part by the National Science 
Foundation under DMR-0130395 and DMR-0325218 and the U.S. Department of Energy 
under DE-FG02-89ER45391 and DE-FG02-03ER15491. The calculations were performed 
at the Minnesota Supercomputing Institute and at the National Energy Research Scientific 
Computing Center (NERSC). 
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TABLE I: Excitation energies of the lowest-energy neutral excitations in benzene. The BSE exci- 
tation energies were obtained with a static kernel. Energies in eV. 





TDLDA 


BSE 


Exp. [24] 


triplet Bf u 


4.5 


3.6 


3.9 


singlet B\ u 


5.4 


4.9 


5.0 


B\ u 


6.2 


6.1 


6.2 


E lu 


6.9-7.2 


7.2 


6.9 



TABLE II: Excitation energies of the lowest-energy neutral excitations in azobenzene. Energies in 
eV. 







TDLDA 


BSE 


Exp. 


TAB 


nir* 


2.1 


2.5 


2.80 [25] 




TTTT* 


3.3 


4.0 


3.94 [25] 


CAB 


nir* 


2.2 


2.6 


2.86 [26] 




TTTT* 


3.2 


3.6 


4.38 [26] 
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Figures 

Fig. ^ Absorption spectrum of benzene, calculated within TDLDA (upper panel) and BSE 
(lower panel). A Gaussian broadening of 0.15 eV was used below 10.0 eV, and 0.5 eV 

n 

above that frequency. The measured spectrum |17| is shown in dashed lines. 

Fig. El Potential energy of the isolated azobenzene molecule in its ground state and the 
two lowest spin-singlet states, obtained within TDLDA (upper panel), and BSE (lower 
panel). The molecular structure is schematically depicted at the two stable configu- 
rations (TAB and CAB) and one intermediate configuration along the isomerization 
path. 
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